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Abstract 

A generic valon potential invariant under U(l) x 50(3) is used to determine whether the hor- 
izontal symmetry of leptons is A4, S4, or some other finite subgroups of 5*0(3). Valons in the 
potential are assigned 50(3) spins 0, 1, and 2, for these are the only ones that can couple to the 
three generation of fermions with horizontal spin 1. This potential causes a breakdown into three 
possible phases with three different symmetries. Phase I has an A4 symmetry suitable for leptonic 
mixing. Phase II has an 50(2) symmetry and phase III has a Z2 x Z2 symmetry, both capable of 
describing Cabibbo mixing of quarks. Phase II has to be rejected on phenomenological ground, but 
phase III yields block diagonal and hierarchical mass matrices appropriate for quarks. No other 
non-abelian symmetry such as 54 is present. 



I. INTRODUCTION 



Symmetry is the foundation of particle physics, so a new symmetry might provide the 
right key to a door beyond the Standard Model. One candidate is a symmetry of the 
three fermion generations which interact identically under the Standard Model. Such a 
(horizontal) symmetry is much discussed, but its existence is not confirmed. Owing to the 
large mass difference and mixing, if it exists it has to be badly broken. That makes it difficult 
to be certain of the existence of a symmetry, much less what it is. While mass difference 
makes it difficult to locate the symmetry, amusingly mixing pattern does help to find it [1]. 

The hunt for horizontal symmetry received an important stimulus when neutrino mixing 
was found to be so regular that it could be described by a tri-bimaximal mixing matrix 
(TBM) [2]. It was also found that this pattern of mixing could be explained by models 
with an A4 horizontal symmetry [3] . A4 is the smallest non-Abelian group that has a three- 
dimensional irreducible representation, making it attractive because most of the existing 
symmetries are in some sense described by the smallest possible groups. 

Nevertheless, the natural symmetry associated with TBM is the group S4 [1], not its 
subgroup A4. The reason why A4 works in those models can be traced back to a dynamical 
choice, allowing an accidental symmetry to promote A4 to S4. From the symmetry point 
of view, S4 is a more natural group, but it is not the smallest group that A4 is. Is the 
leptonic symmetry A4, or S4, or something else? One possible way to decide is to appeal to 
dynamics. 

Both A4 and S4 are finite subgroups of 5*0(3), so if the original unbroken horizontal 
symmetry is S0(3), dynamics of its Higgs valons might determine whether SO (3) breaks 
down to At, or £4, or some other finite subgroup, thereby telling us what the preferred 
symmetry of leptons is. It is this possibility that we want to investigate in the present 
article. 

We shall refer to the quantum number of this horizontal SO (3) as horizontal spin, but the 
word 'horizontal' is often omitted for simplicity because we do not discuss any other kind 
of spin or symmetry here. With three generations of leptons, it is natural to assign them a 
horizontal spin 1. Valons that couple to them in Yukawa terms must then have horizontal 
spins 0, 1, and 2. If SO (3) indeed breaks down to S4 or A4 in the lepton sector, with the 
proper vacuum alignments neutrino mixing is automatically TBM, then the three Yukawa 
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coupling constants can be used to fit the three fermion masses. 

Unfortunately it is known that whatever the dynamics is, it is impossible for SO (3) to be 
spontaneously broken down to either A A or S4, as long as only valons of horizontal spins 0, 
1, and 2 are involved [4]. The reason will be reviewed in more detail in the next section, but 
essentially it is because valons of spins 1 and 2 cannot contain a state invariant under either 
A4 or S4. However, in the case of A4, there is a spin-2 valon state which differs from an 
invariant state under an A A transformation only by a phase factor. No such state exists for 
S4. It is therefore conceivable that in a theory invariant under U (1) x 5*0(3), the extra U(l) 
may provide a compensating phase factor to render this state invariant under A 4 , though 
such an invariant state is still not present in S4. This is indeed the case as we will show 
later, so in this way the dynamics of U(l) x SO (3) chooses A4 to be the appropriate leptonic 
symmetry group over S4. The basic reason for that to be so is because A4 is a smaller group, 
so a state invariant under the subgroup A A may not be invariant under the full group S4. 

We shall refer to the quantum number of this U{1) as horizontal charge. 

The extra symmetry U{1) can be understood as providing a second protection for fermion 
masses. Isospin and hypercharge provides one protection because the left-handed and right- 
handed leptons have different hypercharge and isospin, preventing a mass term to be formed. 
Similarly, the valons here possess a non-trivial horizontal charge, forcing the horizontal 
charges of the left-handed and the right-handed fermions to be different, thereby also for- 
bidding a mass term to be formed. In this way a second protection of the fermion masses is 
obtained. 

A serious impediment to the establishment of a credible horizontal symmetry is the 
qualitative difference between quark and lepton mixings. Instead of having two large mixing 
angles as is the case for neutrinos, quark mixings are all small. The regularity of quark 
mixing is not given by anything resembling a TBM pattern, but merely by a hierarchical 
structure. Quark mixing is dominated by the mixing of the first two generations, with 
the mixing of the second and the third generations quite a bit smaller, and that between 
the first and the third generations the smallest. In the approximation of neglecting the 
smaller mixings with the third generation, the mixing matrix is block diagonal, with the 
third generation unmixed. This means that the mass matrices of the up- and down-quarks 
can be simultaneously put into a block-diagonal form diag(m 12 , m 3 ), where m 3 is the mass of 
the third generation, and m 12 is an arbitrary 2x2 matrix from which Cabibbo mixing of the 
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light masses are determined. We shall refer to a block-diagonal mass matrix of this form as 
hierarchical if |m 3 | is much larger than the matrix elements of mi 2 . From the mass spectra of 
quarks, the mass matrices of the quarks seem to be not only block-diagonal approximately, 
but also hierarchical. 

If horizontal symmetry is really fundamental, how come lepton mixing is governed by a 
symmetry such as At, but quark mixing relies only on a hierarchical structure without any 
apparent symmetry? Is it possible that this qualitative difference is just an illusion at low 
energy, with both of them descending from a common symmetry at high energy? 

Before proceeding along this line it should be mentioned that models can be constructed 
in which both quark and lepton mixings are based on At, but at the expense of giving up 
exact TBM. We will not pursue that possibility here. 

Interestingly the group U(l) x 5*0(3) used to decide the leptonic symmetry to be A 
could also serve as that common symmetry at high energy. This is because the group breaks 
spontaneously into three possible phases with different symmetries, some suitable for leptons, 
and others suitable for quarks. Phase I has a tetrahedral symmetry A that can be employed 
to explain leptonic mixing, and Phases II and III can produce hierarchical block-diagonal 
mass matrices that are appropriate for quark mixing. However, Phase II fails because the 
mixing is between a heavy and a light quark, but Phase III is perfectly viable. 

The details of the U(l) x SO (3) dynamics and the symmetry of the three phases will 
be dealt with in Sec. III. The fermion mass and mixing matrices in different phases will 
be discussed in Sec. IV. There are also two appendices containing mathematical details too 
cumbersome to be discussed in the main text, and a concluding section V. 

II. HORIZONTAL SPIN 

Suppose SO (3) is a horizontal symmetry at some high energy scale. It was mentioned in 
the Introduction that no valon potential can break it down spontaneously to A or S4 if the 
valons have a horizontal spin less than 3 [4] . The reason for that is reviewed in this section. 

If SO (3) can be spontaneously broken down to a subgroup Q, then the valon must contain 
a state invariant under Q. Using the known characters of SO (3) and Q, we can work out 
what irreducible representations (IR) of Q that a horizontal spin- J valon contains. The 
result is shown in Table 1 for J = 0, 1, 2, 3, and Q = A, S4, A 5 . 
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Q 


A 4 


5 4 




IR 


1 V 1" 3 


1 V 2 3' 3 


1 4 5 3' 3 


J = 


10 


1 


1 


J = 1 


1 


1 


1 


J = 2 


111 


110 


10 


J = 3 


10 2 


10 11 


10 10 



Table 1. Multiplicity of horizontal spin- J state in irreducible representations (IR) of Q 
Let us see what the table tells us. 

1. Other than the trivial J = state, the only state invariant under Q (transforms like 
1) is the J = 3 state, with Q = A 4 . There are no invariant states for J = 1 or J = 2. 

2. Although there is no invariant state for J = 2 and Q = A 4 , two of the states (1', 1") are 
pseudo-invariant, in that under an A 4 transformation, these states only gain a phase. 
There is no state in Q = 5 4 nor Q = A 5 which has that property for J < 2. 

3. This observation gives us hope that if we enlarge 5*0(3) to U(l) x 50(3), an additional 
phase from U (1) can be incorporated to cancel the phase of the pseudo-invariant states, 
thereby turning them into invariant states of A4. If that is indeed the case, then A4 
could be a descendent of U(l) x 50(3), though not 5 4 , nor A 5 . 

III. U(l) x SO(3) HIGGS DYNAMICS 

It will be shown in this section that a generic U(l) x 50(3) potential of valons with 
horizontal spins and 2 has three solutions. These solutions differ from one another in 
order parameters and in symmetry. 

Let = (a 2 , a±, a , a_i, a-2) be a complex valon field with horizontal spin 2, and b a 
complex valon field with horizontal spin 0. The subscript of a indicates the horizontal 
angular momentum along some fixed z-axis to be defined later. Let tt(£) be a member of the 
U(l) group obeying w(£i)w(£ 2 ) = +£2), and let w(£)0 = u(£)(f>* = e _ ^0*, u(£)b = 
e^b, u(^)b* = e~ l ^b*. In other words, in additional to the horizontal spin, we also assign cf> 
and b a 'horizontal charge' +1 and 0* and b* a horizontal charge —1. 
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A real renormalizable Higgs potential of 0,0*, 6, and 6* invariant under U(l) x 5*0(3) 
may contain only 0*0, b*b, 0*0*00, 6*6*66, 6*6*00, 0*0*66 types of terms. If we choose 
to couple first 00, and separately 0*0*, before putting them together in the quartic terms, 
then the most general Higgs potential can be seen to be 

V = V (P + V b + Vp,, 

V<t> = (<P<P\9oVo + 92V 2 + 34^4100) - /i 2 (0|0), 

V b = p(b*b) 2 -u 2 (b*b), 

= r{<P\<j ) )b*b + a*b* 2 e + ab 2 e\ (1) 

where Vj is the operator projecting into states with horizontal spin J. The reality of V 
implies that the coupling constants gj,p,r,fi 2 , v 2 are real. The quantity 

6 = £ (-) m a m a_ m (2) 

m=-2 

is proportional to the component of the compound state |00) with zero horizontal angular 
momentum. See (7) below. The potential V^, without the p 2 term, has been used to study 
o?-wave superconductivity, and Bose-Einstein condensates of spin-2 cold atoms [5]. 

V<j, in (1) can be expressed in a more explicit form by getting rid of the projection oper- 

— * — * 

ators. If J\ and J 2 are the horizontal angular momentum operators of valons 1 and 2, then 
one relation between the projection operators can be derived by considering 

(#|(Ji + J 2 ) 2 |00) = (#|07>o + 6P2 + 2OP 4 |00) = 12 + 2(017x10)! • (0| J 2 |0) 2 . (3) 

— * — * 

Choosing the horizontal-spin condensate (0|J|0) along the z-axis, we can write (0| Ji|0)i • 
(0|J2|0)2 = (0|^z|0) 2 '■= (Jz) 2 - Together with the completeness relation Vo + V 2 + V4 = 1, 
we can solve V 2 and V4 in terms of V to get 

^4 = ^(3 + (J 2 ) 2 + 3P ), V 2 = ^(4-(J z ) 2 -10Vo). (4) 
Substituting (4) into (1), and writing 5 (00 (Pol 00) = |#| 2 , the potential becomes 

V 4> = aK 2 + (3{J Z ) 2 + 7 |#| 2 - fi 2 K, (5) 

where 

1 1 12 

« = ^(3^4 + 4^ 2 ), 0=^(94-92), 7 = ^(#0 -9a) + ^(94 -92), (6) 
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and 



^=(0|0> = EM 2 > (^>=E™M 2 > = £(-) ro a_ m «m- (7) 

TO 

No attempt has been made to include couplings with a spin-1 valon in the potential to 
keep it within manageable complexity. We will simply assume that the presence of a spin-1 
valon does not affect the main conclusions reached without it. This would be the case for 
example if the coupling with it is weak, or that its natural energy scale is much lower than 
that of the spin-2 or the spin-0 valons. Moreover, as we will see later, the phase diagram 
obtained from spin-2 valons is not altered by its coupling with a spin-0 valon, so maybe the 
same thing is true for a coupling with a spin-1 valon as well. 

A. Equations of Motion 

There are 12 equations of motion, dV/da* m = = dV/da m (m = 2,1,0,-1,-2), and 
dV/db* = = dV/db, consisting of the following six and their complex conjugates: 



a a m + (3{J z )ma m + i(-) m 9a*_ m + \ab\-) m a^ m = , (8) 

a := a- fi 2 /2K + r\b\ 2 /2K ; (9) 

p \b\ 2 b + a*9b* = , (10) 

p :=p-v 2 /2\b\ 2 + rK/2\b\ 2 . (11) 



The cubic equations in a m , a m , b and b* are more easily solved by first turning them into 
quadratic equations of the order parameters K, (J z ), 9, 9*, and \b\ 2 . This can be achieved by 
multiplying (8) by a m , and separately by (— ) m a_ m , then summing over all m. In this way 
we get 

a K 2 + P(J z ) 2 + 1 \9\ 2 + lab 2 9* = 0, and (12) 

[(a + 7)#+^ 2 ]^ = 0. (13) 

It follows that ab 2 9* must be real. 

Similarly, by multiplying (10) by b*, we get 

p \b\ 4 + a*9b* 2 = 0. (14) 
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Substituting (12) and (14) into (5) and (1), a simplified expression for the Higgs energy 
emerges: 

V = -\{p 2 K + u 2 \b\ 2 ) + ](ab 2 9* + a*b* 2 9). (15) 

B. Solutions 

Solutions can be classified by the order parameters {J Z ),9,K, and \b\ 2 . First, assume 
b^O. 



b^O 



Consider separately the cases 9^0 and = 0. If 9 = 0, then (14) implies p = and 
(13) implies a = 0. Using (11), we obtain 

2p\b\ 2 = v 2 - tK. (16) 

Now consider one after another the two cases, (J z ) = and (J z ) ^ 0. 

I. 9 = 0, (J z ) = 0. From (8), a = 0, hence from (9), 

2aK = p 2 - r\b\ 2 . (17) 

(16) and (17) can be used to solve \b\ 2 and K to get 

+ K= 2jH^. = Ki (lg) 

— r 2 + 4ap — r 2 + 4ap 

Higgs energy can then be obtained from (15) to be 

v = Vi:= _ P» 4 + ^-^\ (19) 
-r 2 + Aap 

One solution for that gives rise to 9 = (J z ) = is 

^ = ^(1,0,0,^,0). (20) 

It will be shown in Appendix A that this solution is unique up to a (7(1) x 5*0(3) 
transformation. 
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II. 9 = 0, (J z ) 7^ 0. (8) implies one and only one a m is non-zero, and that this cannot 
be a for otherwise (J z ) = 0. It is easy to see that (J z ) = mK, so according to (8) 
and (9), a + (5m 2 = 0, hence 

2(a + f3m 2 )K = p 2 - r\b\ 2 . (21) 

K and \b\ 2 can now be solved from (16) and (21) to give 

m2 - T p 2 + 2{a + m 2 (5)v 2 |L |2 K = 2pp 2 - ru 2 = ^ ^ 



11 -t 2 + 4{a + m 2 (3)p ' 1 111 ' " -r 2 + 4(a + m 2 /3)p 

Higgs energy is then obtained from (15) to be 

= ffl 4 + (a + m^)i/ 4 -T/tV 
7/ ' -t 2 + 4(a + m 2 (3)p ' 1 j 

Note that (17) to (19) can be obtained from (21) to (23) simply by setting m — 0. 
The alignment of the spin-2 valon is 

hi = yfK{l, 0, 0, 0, 0), <j>' n = y/K(0, 1, 0, 0, 0), (24) 

for m = 2 and m — 1 respectively. The solutions with m = — 2 and m = —1 can be 
obtained from these by a 7r-rotation about the y-axis. 

III. 6 7^ 0. In this case we must have (J 2 ) = 0. Otherwise, as argued before, there is one 
and only one a m ^ 0, and this is not a . As a result, according to (2) or (7), we must 
have 9 = 0, contradicting the assumption that 9 ^ 0. 

The solution in this case is considerably more complicated because we are no longer 
forced to have a — 0. 

Multiply (13) by 9* , then compare it with (12) and with (14), we get 

a,K 2 + 1 \9\ 2 = poN 4 A (25) 
(ao + 7 )|0| 2 = Po\b\ 4 /2, (26) 

which shows that \9\ 2 = K 2 . It also shows that p and a + 7 are both positive, or 
both negative. 

Let us first assume both of them to be negative. Then taking the absolute value of 
(13) and (10), and use (9) and (11), we get 

(a, + 1 )K = (a + 1 )K+ l -(r\b\ 2 -p 2 ) = -\\a\\b\ 2 , (27) 

Po \b\ 2 = p \b\ 2 + 1 (rK -v 2 ) = -\a\K. (28) 
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These two yield the solution \b\ 2 = \bni\ 2 and K = K in , with 

-(r+ \a\)u 2 + 2pp 2 



2(a + 7 )z/ 2 -(2| ( x 




4(a + 7)p- (r + 


a\)(r + 2 





\hn\ 2 

They cause the Higgs energy in (15) to be 

V : = Vm 
iVi 



K 



in 



4(a + 7)p- (r+ |<r|)(r + 2 



p 4 Ni + /i 2 z/ 2 A" 2 + v 4 N 3 
[4(a + 7)p- (r + |<r|)(r + 2|a| 
p[4p(a + 7 ) - r (r + 2H)], 



n 2 ' 



N 2 = (r + 2\a\) -Ap{a + 7) + (r + 



-(r + |(x|) 2 (a + 7) + 4p(a + 7) 5 



(30) 



If both p and a + 7 are positive, then the right-hand side of equations (27) and (28) 
change a sign. In that case, solutions (29) and (30) are still valid provided we change 
\a\ in every expression to —\cr\. 

One spin-2 vacuum alignment that gives rise to \9\ = K and (J z ) = is 

4>ni = VK(smi]/y/2, 0, cost], 0, smrj/y/2), (31) 

where 77 is an arbitrary angle between and ir. It will be shown in Appendix A that 
all other solutions can be obtained from this one by a U(l) x 5*0(3) transformation. 

For convenience, these three type of solutions with b 7^ are summarized in Table 2. The 
last column indicates the condition under which the solution is valid, and the numbers within 
parentheses are equation numbers describing the solutions. 



sol 


\e\ 


(J z ) 


\b\ 2 


K 


V 





cond 


I 








(18) 


(18) 


(19) 


(20) 


a = 


II 





mK 


(22) 


(22) 


(23) 


(24) 


(7 = 


III 


K 





(29) 


(29) 


(30) 


(31) 





Table 2. A summary of the b 7^ solutions 



b = 



It is clear from (10) that there are solutions with 6 = 0. It is also clear from (8) that the 
solutions for spin-2 valons are identical to those obtained above, after setting a = r = 0. 
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From (15), the potential V is also identical to those obtained before by setting a = r = 
u 2 = 0. 



C. Energy and Phases 

Stability demands a condensate to have the lowest energy. It is therefore necessary to 
determine which of the solutions in Table 2 fits that role. The result turns out to depend 
only on the relative sizes of the parameters a, (3, and 7, as shown in the phase diagram in 
Fig. 1, and not on the parameter r. Here is how that conclusion is arrived at. 

First of all, if a ^ 0, then the only solution and the only phase is III, so there is no phase 
diagram and no need for a further discussion. From now on, we will assume a = 0. 

The parameters /i and v have the dimension of energy, and the rest are dimensionless. 
In the absence of couplings between the (horizontal) spin-2 and spin-0 valons (o = r = 0), 
they define the energy scales of the spin-2 condensate K and the spin-0 condensate |6| 2 , 
respectively. The spin-0 condensate is given by \b\ 2 = v 2 /2p, and the spin-2 condensate is 
given by either I. K = p 2 /2a, II. K = p 2 /2(a + m 2 /3), or III. K = p 2 /2(a + 7), depending 
on which of the solution has the least energy. Note that since \b\ 2 and K are non-negative, 
the parameters p, a, (a + 7), (a + 4/3) must also be non-negative. 

The energy V in this un-coupled scenario is given by (15) to be V = (p 2 K + u 2 \b\ 2 )/2, 
namely, 



depending on which of the energies is the smallest. If /3, 7 > 0, the lowest energy occurs in 
phase I. If < 0, II has a lower energy than I, and m = 2 always has a lower energy than 
m=l, hence m — 1 is never the condensate and we will consider it no further. If 7 < 0, 
then II is still the condensate if A\f3\ > I7I, otherwise III is the condensate. The phase 
diagram based on this discussion is shown in Fig. 1. 

If r is not zero, the formulas are complicated and the detail is relegated to Appendix 
B. Nevertheless, it turns out that Fig. 1 remains valid for all r, at least for those solutions 




(m = 1,2), or 



(32) 
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a+y=0 



p 



III 



I 




7 



a+4fi=0 



FIG. 1: Phase diagram 



defined and connected continuously to the un-coupled solutions. 
D. Symmetry of Solutions 

Each of the three solutions in Table 2 is invariant under some symmetry group Q C 
U (1) x 50(3). One way to determine the symmetry is to use Wigner's D-matrix to compute 
Y, m ' a m'T^'!r7mi. a ' i P' i an d equate it to a m e~ lp ' for some U (1) phase e~ ip ' . The Euler angles 
of rotations are denoted here by a', ft, 7', with a prime to tell them apart from the coupling 
constants in (1) and (6). This gives us five equations (—2 < m < 2) to solve for four 
unknowns a!,ft,^',p'. If there is no solution other than a' — ft — 7' = 0, then there is no 
non-trivial symmetry. Otherwise, the solution of (ot',ft,^') yields the 5*0(3) element and 
e ip ' the U{1) element in Q. 

I. Tetrahedral Symmetry 

The solution (p I = ^K/3(l, 0, 0, y/2, 0) has a tetrahedral symmetry. It is invariant under 
the two i7(l) x SO (3) operations 



-ifioJy -iirJ z ~ ini/3 -2mJ z /3 



(33) 




Their explicit spin-2 matrix representations 
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are gi = u 2 V J ~ 2 (0, 0, 2tc/3) = diag(l, u, u 2 , 1, ou), and 

/ 1 2^2 2^/6 4^2 4 \ 
2^2 5 2\/3 -2 -4^/2 



g I = V J=2 (O,0 o ,n) = 
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2\/6 2^ -3 -2^/3 ^6 
4^/2 -2 -2^/3 5 -2^/2 
V 4 -4^/2 2^/6 -2>/2 1 / 

They can be used to verify the invariance gi<pi = <pi and gi<pi = <pi- 
These two operators satisfy 

g 2 i = fi = {gigif = i, 



(34) 



(35) 



so they generate an A A group, the invariant group of a tetrahedron. Note that the U(l) 
phase factor in (33) is not necessary for gj to be a generator of Qi = A4, but it is necessary 
for (pj to be invariant under gi. 

II. Ferromagnetic Symmetry (m = 2) 



The solution <\> u = vK(l, 0, 0, 0, 0) of (24) is invariant under the operation g n = 
u(2£)e-'^ h for any £, hence Q n = SO(2). 

Originally there were two additive quantum numbers, the horizontal charge Q and the 
third component of the horizontal spin m. In the ferromagnetic phase, only m! = m — 2Q is 
conserved. This is reminiscent of the Standard Model, in which both hypercharge and weak 
isospin are broken in the Higgs phase, leaving the electric charge to be the only additively 
conserved quantum number. 

III. Nematic Symmetry 

The solution <p in = y/~K(smr}/y/2, 0, cost/, 0, sinrf/y/2) is invariant under g ni = 
e~ mJz , gjn = e~ mJy , and g~ in = e~ mJx . Any two of those three will generate the sym- 
metry group Qui = Z 2 x Z 2 valid for any r\ in (ftni- 

At special 77's, the symmetry group is enlarged to 5*0(2). At rj = 0, 4>m oc (0, 0, 1, 0, 0) 
is invariant under any rotation about the z-axis. At r\ = ir/3, (pm oc (V3, 0, \/2, 0, y/3) is 
invariant under any rotation about the y-axis, and at r\ = 2tt/3, (fini oc (\/3, 0, — y/2, \/3) is 
invariant under any rotation about the x-axis. 
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These symmetries of the three phases are summarized in Table 3. Note that the generators 
need a U(l) factor in land II, but not III. 



1 

sol 


V 





Q 


9,9 


T 

1 




V y (M,0, v2,0J 


A 

A A 




II 




v^(l,0,0,0,0) 


SO{2) 




III 


V 


sin ?7, 0, cos 77, 0, ^ sin 77) 


Z 2 x Z 2 









V^(0,0,1,0,0) 


SO{2) 






tt/3 


(v^3, 0, >/2, 0, V3) 


SO{2) 


g — i£ Jy 




2tt/3 




SO{2) 


e -i£,Jx 



Table 3. Symmetry of the spin-2 solutions 



IV. MASS MATRIX AND MIXING 



We will discuss in this section how to compute mass matrix and fermion mixing for the 
three phases. Phase I with an A 4 symmetry is presumably suitable for leptons. Phases II 
and III with a hierarchical structure are potentially good for the quarks, but phase II turns 
out to be unphysical because its main mixing is between a heavy and a light quark. 

The mass matrix in these phases is most easily computed when valons are expressed in a 
matrix form. The following subsection shows how this can be done for valons of horizontal 
spins 0, 1, and 2. 



A. Valons in matrix form 



In the Cartesian basis, a spin-2 valon is described by a traceless symmetric matrix (or 
tensor) $ = = 1, 2, 3 or x,y, z), whose entries are given by the components of the 

spin-2 state = (a 2 , a±, a , a_i, a_ 2 ) to be 

/ |(a 2 + a_ 2 ) - ^a §(a_ 2 - a 2 ) |(a_i - a x ) \ 

;(a_ 2 -a 2 ) -\{a 2 + a_ 2 ) - ~^a §(a_i + ai) 

Ua-i + ai) 



V 



a_i — ai 



a 



(36) 



The normalization is chosen so that K := J2 m |^m| 2 = Tr($^$). With that normalization, 
the scalar 9 is given by 9 = I] m (— ) m a m a_ m = Tr($ 2 ). The three components of the vector 
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(J) is proportional to the three independent components of the anti-symmetric hermitian 
matrix with (J 2 ) = (J z ) 2 = (E m m\a m \ 2 ) 2 = -2Tr([$t, $] 2 ). 

Valons of spin and spin 1 may also be expressed as matrices. The matrix form B of a 
spin-0 valon b, and T of a spin-1 valon c = (ci, c 2 , c 3 ) in the Cartesian basis, are 



B = 



(b 





°^ 









C3 


-C2\ 





b 









-C3 





Cl 


\o 











C2 


-ci 


o / 



(37) 



2 .. 



Normalization is chosen so that Tr(B^B) = \b\ 2 and Tr(r^r) = J2i |q 

Under an 5(9(3) transformation, A — >■ ft^A/i for A = $, £?, and T, where ft, is the three- 
dimensional representation of the 5*0(3) rotation in the Cartesian basis, obtained from 
T> J=1 (a', /3',7') through a similarity transformation that converts spherical basis into Carte- 
sian basis. 

In particular, the matrix expressions of the three condensates (pi, <pu, and (pm are 



1 K 
2V 3" 



ii 



/ 1 -i V2\ 

-i -1 a/2z , $ 
\\/2 \/2i / 
/ — cos 77 + \/3 sin r\ 

in \j — cos// \'.'5sin// . (38) 

V o 

For a generic rj, §ui is invariant under the rotations generated by any pair of the inver- 
sions x — > —x, y — > —y, z — > —z, reflecting its nematic structure Z 2 x Z 2 . 

For and § n , the explicit form of hj,hj, hu, the three-dimensional Cartesian repre- 
sentations of the S0(3) operators gi = e 
are 



/ 1 -i 0\ 

-i -1 

V o oo/ 

\ 

cos rj — \/3 sin 77 

2 cos n / 



/ 1 2 v / 2\ 



hi 



hi 



0-3 

\2 V / 2 -1 y 

Using them, the correct transformations 



f-1 -V3 0\ 

v/3 -1 
V 2/ 



/ cos £ — sin £ \ 
sin £ cos £ 

V 01/ 



(39) 



h\$ I h I = & I , h\$ jhj = u$ 7 ; h\j® n h H = e 2l? $jj 



(40) 



can be verified to be true. 
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B. Fermion Mass Matrices 

With the left-handed and right-handed fermions taken to be horizontal vectors, the mass 
matrix M is just the expectation values of valons they couple to, expressed in the matrix 
form: 

M = $ + 5 + r. (41) 

The fields a m ,b,Ci in this formula should be interpreted as expectation values. We will 
assume they can be assigned any value that suits the phenomenology. No Yukawa coupling 
constant appears because they have been absorbed into the fields. 

In order to gain a qualitative understanding of the three phases, let us examine the 
simplest case when the expectation values are given only by the condensates in the last 
section. Namely, $ is given by (38), and B = T = 0. 

The squared-mass spectra in these three phases can then be obtained from the eigenvalues 
of M := MtM = $t$ to be 

t K K K 

3 ' 3 ' 3 • 

II. 0, 0, K 

III. K(s — c) 2 , K(s + c) 2 , 4Kc 2 , where s = sin rj/ a/2 and c = cos r)/\/6. 

On the surface solution II seems to be tailored made for quark mixing: the mass matrix 
(38) has a block diagonal form and the mass spectrum is hierarchical, with one heavy mass 
K and two light masses 0. The trouble is, one of the light masses comes from the unmixed 
33 entry. Since both the up- and the down-quark mass matrices have these properties, quark 
mixing occurs only in the 1-2 block, giving rise to a mixing between a heavy quark and a light 
quark, not between two light quarks. The inclusion of B and Y will not change the behavior 
qualitatively, so phase II has to be rejected on phenomenological grounds. However, the 
remaining two phases will be put to good use as we shall see in the next two subsections. 
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C. Quark Sector (Phase III) 



In order to force the mixing to be between the first two generations, we will assume only 
c 3 7^ when the vector acquires an expectation value. In that case the mass matrix is 

/y/K(-s + c)+V c' 3 \ 

M = $ 7// + B + T = -c 3 -VK(s + c) + b' , (42) 

V 2^/Kc + b 1 ) 

where c = cos^/v^, s = sin^/y^, b' = b/y/3, and c 3 = c 3 /\/2. This matrix has the correct 
block-diagonal form. It can be made hierarchical by adjusting the three parameters K,b,r] 
so that the 33-entry is much bigger than the 11- and 22-entries. We will assume 77 to be in 
the second quadrant so that the hierarchical mass spectrum mi < m 2 < m 3 can be obtained 
with b' and s positive, c negative, and \s\ ^> |c|. We shall also assume c' 3 := 27 to be purely 
imaginary. In that case M is hermitian and it can be diagonalized by a unitary matrix. 

Since the third generation stands by itself, we need to consider only the submatrix m i2 
in the upper-left block. It can be diagonalized using the formula 

(cos 9 —ism9\/ a ij\ ( cos 9 ism9\ ( m\ 0\ 
— i sin 9 cos 9 ) \ —ij f3 ) \ i sin 9 cos 9 ) \ m 2 / 

where 

mi + m 2 = a + (3 = 2{b' - y/Ks), m 2 -m x = J(/3 - a) 2 + A-f 2 = . 

tan2^ = ^ , (44) 

The constants K,b,r) can be fixed from the three quark masses mi,m 2 ,m 3 . Adjusting the 
relative size between 7 and (3 — a, the angle 6 can be chosen to be anything we want. The 
diagonalization matrix is unitary but not orthogonal, but it differs from an orthogonal matrix 
just by some column and row phase factors which we are allowed to insert by adjusting the 
phases of the up and down quarks. Thus the Cabibbo angle is just the difference of the 
^-angles for the down and the up quarks, something that is up to us to adjust. 

This simple quark model is presented to illustrate that Cabibbo mixing with hierarchical 
quark mass matrices can be accommodated in the nematic phase III. 
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D. Lepton Sector (Phase I) 



Before plunging into the construction of a horizontal-spin model of leptons in phase I, it 
is useful to recall the representations of A 4 , and how trimaximal as well as tri-bimaximal 
mixings are obtained. 

A4 has four irreducible representations (IR), I, I', 1", and 3. The representation of its 
generators gi and gi, satisfying the relations gj — gj — (gigi) 3 = 1, is given in Table 4 in 
the basis where gi is diagonal (u> = e 27 ™/ 3 ). 



IR 


1 


V 


1" 


3 


9i 


1 




CO 2 


diag(l, u, u 2 ) 












9i 


1 


1 


1 















Table 4. Irreducible representations of A A 



u : = 



1 

71 



We shall refer to the basis used in Table 4 as the l F basis', and denote the 3-representation 
of gi and gi there as ji and tu 2 ji. They differ from the corresponding SO (3) elements hj 
and hj in (39) only by a unitary transformation, 

/l i \ 

1 —i , hi — u^jju, hi = v)jiu. (45) 
\0 V2J 

When both the left-handed and the right-handed fermions are A 4 triplets, the Yukawa 
valons coupled to them belong to 1, 1', 1", 3 S , or 3 a , where 3 S is a triplet that couples sym- 
metrically to the two fermions, and 3 a a triplet that couples anti-symmetrically. To obtain 
trimaximal mixing, we should assign the expectation value of a valon in the charged-lepton 
(neutrino) sector to be an invariant eigenvector of gi (gi), for every irreducible representa- 
tion. If that is not possible then the expectation value shoud be taken to be zero. With this 
assignment trimaximal mixing is guaranteed. Once a trimaximal mixing is thus obtained, 
an appropriate adjustment of the Yukawa coupling constants can promote the mixing into 
a tri-bimaximal mixing. 

To implement this procedure in the U (1) x 5*0(3) theory where gi = e~ tl3 ° Jy e~ lnJz and gi = 
u 2 e~ 2mJz ^ 3 , it is necessary to identify the nine A 4 states in representations 1, 1', 1", 3 S , 3 a 
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from the nine horizontal states endowed with horizontal spins 0, 1, and 2. The spin-0 and 
spin-1 valons are B and Y in (37). As to the spin-2 state, we assume it either comes from 
the dynamical valon $':=$ — $/ acquiring a vacuum expectation value at a lower energy, 
or a new low-energy spin-2 valon which does not participate in the high energy dynamics 
that determines the phase structure. To simplify writing in the first scenario, which we will 
implicitly assume from now on, the prime will be omitted. Instead, we will replace $ in (41) 
by + $. Equivalently, in the spherical basis, is to be replaced by 0/ + 0. 

It is obvious from Table 4 and the explicit form of gi,gi that the spin-0 state is a 1" 
of A 4 . The spin-2 state can be decomposed into a mixture of 1,1', 3 of A 4 as follows. 
We already know that (1, 0, 0, v^2, 0) transforms like 1. Using (34) it is easy to verify that 
(0, — y/2, 0, 0, 1) transforms like V. After projecting these two out, what remains transforms 
like 3. In this way the decomposition of a spin-2 state = (02, a±, a , a_i, a-2) is obtained. 

0i = -(a 2 + V2a_i) (1,0,0,^,0) := ^(1,0,0,^,0) , 

0i' = ^(a_ 2 - v / 2ai) (0,-^,0, 0,l) := VZd (0, -y/2, 0, 0, l) , 



2a 9 



_i, ai + V2a-2, 3a , a_i — V^a^ 2a_ 2 + V^- 



a,i 



= 01 + 01' + 03- 



(46) 



03 is a 5-dimensional vector with 3 independent components. It can be projected into the 3 
space to get 

1 



v = (vi,v 2 ,v s ) = -W-ai - v^a-2, v^ao, + ^202) := 0"P- 

v3 



(47) 



The components of t> are chosen to commute with the 5*0(3) rotations, so that (4>gi)V = 
{4>V)fi and (4>gi)V = (<f)V)u 2 fj, where // and // are the SO (3) part of the spin-1 A 4 
generators in the spherical base, namely, // = diag(o; _1 , l,w) and 

l\ 2 2 \ 

fl = ~\ 2 1 2 ■ ( 48 ) 
\2 2 1 J 

These two differ from the 3 representation ji = gi and ji = ^ 2 gi = diag(cu, cu 2 , 1) in the 
last column of Table 4 only by an inconsequential unitary transformation W, 

( \\ 



W := 



-10 
\ 1 0/ 



Ji = WfjW\ jj^WfjWl 



(49) 
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c 



u 



1 

71 



They also differ from hi and hi of (39) just by a unitary transformation C that converts 
spherical harmonics into Cartesian coordinates: 

/-It \ 

o s/2 , hi = C^\C, h I = C ] hC. (50) 
\i » o/ 

The normalization of v is chosen so that its norm \v | 2 := X)j K| 2 = |03| 2 := Em |(03)m| 2 - 
Besides the spin-2 valon 0, the mass matrix M also receives contribution from the spin-0 

valon b which is a 1", and the spin-1 valon c which is a 3 a because T is an anti-symmetric 

matrix. 

Next, let us compute the expectation value of valons needed to yield trimaximal mixing. 
Those used in the charged-lepton mass matrix M e are given by the invariant eigenvectors of 
gi, and those used in the neutrino mass matrices M v are given by the invariant eigenvectors 
of gi. In case the operator does not have an eigenvalue +1, the corresponding expectation 
value is zero. The solution to this requirement is summarized in Table 5. Normalization 
cannot be determined so those listed in Table 5 are really vacuum alignments, computed in 
the F-basis where gi is diagonal. The vectors vf and cf are vectors v and c expressed in 
that basis. The table also gives relations between a m and between q needed to obtain such 
vacuum alignments. They can be computed either from the invariant eigenvectors vp and 
Cp, or more directly from the condition that ~h\~M~h x = uM e and h\M v h x = M v . 



valon 


e 


d 


vf 


b 


c F 


IR 


1 


V 


3 


1" 


3 




1 





(1,0,0) 





(1,0,0) 








CLi = do = d-2 = 




Cl = c 3 = 


M v 


1 


1 


(1,1,1) 


1 


(1,1,1) 








— \f2a 2 + ai + a_i + = 
— v^2a 2 + \/3a + a_i = 




ci = V2c 3 
c 2 = 



Table 5. Vacuum alignment of horizontal valons 

Using Table 5, we obtain the F-basis mass matrix «($/ + $ + B + T)v) to be 

/ a' 2 \ 

a , _ 1 /V2-V2c 2 

\a'_ 1 /V2 + >/2c2 / 



M e = 



20 



/ -v' 3 + b' - ic 3 e' + 2v' 3 d - v' 3 + ic 3 \ 
M v = d + 2v' 3 b'-v' 3 + ic 3 e'-v' 3 -ic 3 , (51) 

\ e' — v 3 + ic 3 d — v' 3 — ic 3 b' + 2v' 3 J 

where k = ^Jk/3, a! 2 = a 2 + n, a!_ x = a_i + \/2k, e' = e + k, b' = b/ y/3, and v' 3 = v 3 / 

In this F-basis, jiM u j\ = M v and j\M e jj = uM e . Although jj is diagonal, M e is not 
because of the phase factor u. However, this phase factor cancels out in M e = M\M e so M e 
is diagonal, with eigenvalues m 2 e = (V^c 2 + a-i/y/2) 2 , m 2 ^ = a|, m 2 = (V2c 2 — a^i/^/2) 2 . 
They can be used to determine the parameters c 2 , a 2 , a_i from the masses of the charged 
leptons. In particular, c 2 and a_i must have opposite signs to yield the correct hierarchy of 
masses. 

The expression for M v in (51) can be used for the Dirac mass matrix of the neutrinos as 
well as the Majorana mass matrix of the heavy right-handed neutrinos, with two different 
sets of parameters of course. This matrix is magic because the sum of every row and every 
column is equal to b' + d + e', hence the mixing matrix is trimaximal [6]. If we choose the 
parameters such that e' — d and C3 = —3iv' 3 , then the mixing matrix is 2-3 symmetric and 
hence bimaximal [6]. With this choice, the neutrino mixing matrix is both trimaximal and 
bimaximal so it has the TBM form. 



V. CONCLUSION 



We propose to resort to dynamics to determine which of the finite non-abelian subgroups 
of SO (3) is the appropriate horizontal symmetry of leptons. We do so by using dynamics 
of valons of horizontal spin not larger than 2, invariant under the smallest group containing 
5*0(3) that can spontaneously break down to at least one non-abelian subgroup. It is found 
that using U (1) x 50(3), the only non-abelian subgroup it can break down to is thereby 
making it the preferred horizontal symmetry of leptons. Other than A4, this dynamics also 
produces two other phases, a phase II with a 50(2) symmetry, and a phase III with a 
Z 2 x Z 2 symmetry, both capable of describing Cabibbo mixing of the quarks with a mass 
hierarchy. It turns out that phase II yields the wrong mass hierarchy and has to be rejected, 
but phase III is perfectly viable. Thus U(l) x 50(3) can also be thought of as a common 
high-energy horizontal symmetry where both the Cabibbo mixing of quarks and the tri- 
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bimaximal mixing of neutrinos originate. Explicit mass matrices in the quark sector (phase 
III) and in the lepton sector (phase I) are constructed to illustrate these features. 

Technical complication prevented us from including a fully coupled spin-1 valon into the 
potential in this paper. CP violation is also left out. To include the latter it is presumably 
necessary to go beyond SO (3) to SU (3). These points are being investigated. There are also 
several other dynamical issues that requires further study. For example, we have assigned 
phase I to the lepton sector and phase III to the quark sector because that seems to be what 
phenomenology demands, but it would be much nicer if a dynamical mechanism can be 
found to force that to happen automatically. This paper is of an exploratory nature, to test 
the idea and to demonstrate the feasibility of a symmetry group common to the quark and 
lepton sectors, a higher group that can even select the lepton symmetry at lower energies. It 
is not meant to present a complete dynamical model so detailed questions such as whether 
this higher symmetry is local or global, and the associated question of Goldstone bosons, 
etc., have all been left out. 

I am grateful to Fei Zhou, We-Fu Chang, and James Bjorken for stimulating discussions. 
Part of this work was carried out while visiting the Theoretical Science Center at Tsinghua 
University in Hsinchu, and the Academia Sinica in Taipei, in the summer of 2010, whose 
hospitality I would also like to acknowledge. 
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Appendix A: Uniqueness of the spin-2 solutions 



The purpose of this appendix is to show that the (J 2 ) = (J z ) 2 = solutions (20) and 
(31) are unique up to £7(1) x 5*0(3) transformations. 

As discussed in Sec. IV, a spin-2 valon can be represented by a traceless symmetric 
matrix $, with (<f>\<j>) = Tr($t$), 9 = Tr($ 2 ), and (J 2 ) = -2Tr([$t, $] 2 ). 

Since [<&' , <3>] is hermitian, $] 2 i s positive semi-definite, with non-negative real eigen- 
values. To have a zero trace needed for (J 2 ) = 0, all eigenvalue must be zero, implying 
[$"!",$] = 0. Therefore $t and $ can be simultaneously diagonalized by an orthogonal 
transformation. 

Suppose diag(a, b, —a — b) is the diagonal form of the traceless matrix $. Applying a 
£7(1) transformation if necessary we may assume a to be real. The 9 parameter is given by 

# = Tr($ 2 ) = 2(a 2 + b 2 + ab). (Al) 

If 9 = 0, then it is still zero if we multiply both sides of (Al) by a or by b. Equating these 
two identities, we get a 3 = —a(b 2 + ab) = —b(ab + a 2 ) = b 3 . Therefore b/a is a third root of 
unity. This root cannot be 1 for then 9 ^ 0, hence it must be to = e 27 ™/ 3 or to 2 = e 47 ™/ 3 . Thus 
the eigenvalues of $ are a, aou, and au 2 . Since a is fixed by the norm K = ((f>\(f>) = Tr(&&) 



to be a = ±y A"/3, the solution $ for 9 = J z = is unique up to a sign and a U(l) x 5*0(3) 
transformation, hence <f>i in (20) is unique up to a 17(1) x 50(3) transformation. 

As a check, we can compute directly the eigenvalues of $/. They turn out to be 
v /a73, \[K/3 oj, and ^K/3 u 2 . 

If \9\ = K, with a C/(l) transformation we can render 9 = K. In that case 9 = a 2 + 
b 2 + (a + b) 2 = K = \a\ 2 + \b\ 2 + \a + 6| 2 . Since a is real, so must be b. As a function 



of b, a reaches a stationary point at a = —2b, whence 9 = 3b 2 = K, giving b = ±yK/3 
and a = =(=2^7^/3. We can now parametrize a to be a = 2^Jk/3cosi]. Substituting this 
into (Al) with 9 = K, we can solve b to get b = ^ K/3(— cost] ± v^sinr?). Thus the 
eigenvalues of $ are \J^K/3 (2 cost], —cosi] + y/3 sin 77, — cost? — v^sin??). Comparing this 
with the expression for $> IH in (31), we see that the 77 here is the same as the 77 there. It 
also shows the uniqueness of ipni up to a £7(1) x 50(3) transformation. 
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Appendix B: Interacting condensates 

The dependence of V, K, \b\ 2 on r and a will be worked out in this appendix for phases 
I, II, and III. We will also work out what the phase diagram Fig. 1 becomes when r ^0. 

For the sake of this discussion, we shall assume the spin-2 condensate in the absence of 
coupling to be larger than the spin-0 condensate, namely, ^Jav 2 j ^fpix 2 := r < 1. 

1. 9 = 0, (J z )=0 

The condensates are given in (18) and the energy in (19). 

First consider how Vi changes with r. At r = 0, Vj = (Vb)j = — (yU 4 /4o; + z/ 4 /4p). 
Since (dVi/dr) T=0 = p 2 z/ 2 /4ap is positive, \Vi\ decreases with increasing r in the vicinity of 
t = 0. Farther away, Vj has two stationary points where dVj/dr vanishes. The smaller one 
is located at r = t_ := 2av 2 /fi 2 := r r, where r := and the larger one is located 

at t + = To/r. Vj reaches a local maximum (V})_ = —fi 4 /Aa at r_ and a local minimum 
(Vj)+ = — z/ 4 /4p at t+. The local maximum is smaller (lower) than the local minimum 
because r < 1. The sum of the values at the local maximum and local minimum happens 
to be the value Vj at r = 0. 

There are two asymptotes located at r = ±r where Vj = ±oo. They divide the graph of 
Vj vs r into three branches. Since < t_ = r r < r , the smaller stationary point appears 
in the central branch. The other stationary point at r + = r /r > r is to be found in the 
right-hand branch. 

The left-hand branch starts at Vj = at r = — oo and rises monotonically to +oo at the 
left asymptote r = — r . The central branch starts at Vj = — oo at the left asymptote, rises 
to a maximum at r = r_, then drops back to — oo at the right asymptote. The right-hand 
branch starts at V} = +oo at the right asymptote, descends to a minimum at r = r+, then 
rises back to at r = +oo. 

The central branch is sketched in Fig. 2. 

Next, let us see how the condensates b 2 and K vary with r. At r = 0, b 2 = u 2 /2p 
and K = /i 2 /2a. The product of these two happens to be the slope of Vj at that point. 
Both of the b 2 - and A'-slopes at r = are negative, b 2 decreases monotonically from the 
value z/ 2 /2p at r = to the value at r = r_, beyond which 6 2 becomes negative and 
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FIG. 2: The center branch of the Higgs potential Vj vs r 

unacceptable. For that reason we need not consider the central branch in Fig. 2 to the right 
of r_. The condensate starts from the value /i 2 /2a at r = 0, decreases to a minimum 
value (fx 2 /2a)(r /2Q at r = r m = (r , where ( : = (1 — Vl — r 2 )/r. From r = r m to r_, 
.fT increases from this minimum back to the value /z 2 /2a. To be consistent, we should have 
C < 1 and r < 2£, which indeed follows from the condition < r < 1. 

For negative r, both b 2 and -ft" increase monotonically to the value +oo at the left asymp- 
tote T = — To- 

The behavior of b 2 , Kj, and Vj is summarized in Table 5. The definition of various 
T a is given in the second row, with r := y/au 2 / \fpfi 2 and £ := (1 — \J\ — r 2 )/r. Also, 
£ := r 2 /2 + r/2( - 1. It follows from < r < 1 that < r/2 < ( < r and < £ < r 2 . 
Hence b 2 and V are monotonic in the range — tq < r < r_, but .fT has a minimum located 
at t = T m . 



T 


-TO 





7*m 


t_ 


def 


-2^/ap 







rr 


b 2 


+oo 


v 2 /2p 


z/ 2 /4p 





K 


+oo 


H 2 /2a 


(/i 2 /2a)(r/2C) 


H 2 /2a 


V 


— oo 


-(/i 4 /4a)(l+r 2 ) 


-(/i 2 /4a)(l+,) 


-fu 4 /4a 
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Table 5. The variation of (18) and (19) with r in the case of 9 = (J z ) = 



2. = 0, (J z ) = 2K 

The condensates in this phase are given in (22) and the energy given in (23), both with 
m — 2. In particular, V u in (23) is the same as Vj in (19) after having a replaced by 
a' = a + 4/3. 

Since 

dV L= {-2pl+j?rf 

da ' (-r 2 + 4ap) 2 " ' 1 ' 

phase II has a higher (lower) energy than phase I if /3 > (< 0). In other words, the 
boundary between phases I and II in the phase diagram of Fig. 1 remains the same even if 
the coupling r ^ is turned on. 



3. |0| = K, (J z ) = 

The condensates in this phase are given in (29) and the energy given in (30). 

If a = 0, then (29) can be obtained from (19) by replacing a with a" = a + 7. On account 
of (Bl), the boundaries of phase III with phases I and II remain the same as Fig. 1 even in 
the presence of r-coupling. 

If a 7^ 0, the r-dependence is very complicated. However, if the spin-0 condensate b 2 is 
non-zero, the solutions in phases I and II do not exist for a 7^ 0, so this phase III remains 
the only phase no matter what r is. 
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